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ABSTRACT 

Context. Observations of the effect of microlensing in gravitationally lensed quasars could potentially be used to study the structure 
of the source on distance scales down to the size of the supermassive black hole powering the quasar activity. 

Aims. We search for the microlensing effect in the y-ray band using the signal from a gravitationally lensed biazar B0218+357. 
Methods. We develop a method of deconvolution of contributions of two images of the source into the y-ray band flaring lightcurve. 
We use this method to study the evolution of the magnification factor ratio between the two images throughout the flaring episodes. We 
interpret the time variability of the ratio as a signature of the microlensing effect and derive constraints on the physical parameters of 
the y-ray source by comparing the observed variability properties of the magnification factor ratio with those derived from numerical 
simulations of the microlensing caustics networks. 

Results. We find that the magnification factor ratio has experienced a change characteristic for a microlensing caustic crossing 
event during a 100 d flaring period in 2012. It has further changed between 2012 and a recent flaring episode in 2014. We use the 
measurement of the maximal magnification and duration of the caustic crossing event to derive an estimate of the projected size of 
the y-ray emission region in B0218+357, Ry ~ lO''' cm. This estimate is compatible with a complementary estimate found from the 
minimal variability time scale. The microlensing / minimal variability time scale measurements of the source size suggest that the 
y-ray emission is produced at the base of the biazar jet, in the direct vicinity of the central supermassive black hole. 
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P 1. Introduction 

^ Almost hundred years after the discovery of the first large scale 
I ^ I jet in an Active Galactic Nucleus (AGN) ( |Curtis|1918 1, the ori¬ 
gin of the jets and the details of the mechanism of their produc¬ 
tion and high-energy activity remain uncertain. 

Radio observations are now starting to reach the angular res- 
olution sufficient to resolve the details of the region of the jet 
formation in the direct vicinity of the supermassive black hole(s) 
in nearby AGN ( [Doeleman et al.|[2012| l. The angular resolution 
of y-ray telescopes will not be comparably high in any foresee¬ 
able future. This precludes the possibility to get the necessary 
imaging data which would provide answers to the most basic 
questions, such as e.g. the location of the site of the y-ray emis- 
sion region in blazars (AGNs with jets aligned along the line of 
sight, |Urry & Padovani|1995|). 

The only chance to overcome the limitation stemming from 
the limited angular resolution of y-ray telescopes is provided 
by the effect of the gravitational microlensing observable in 
the blazars / AGN which are parts of the strong gravitational 
lens systems ( [Chang & RefsdaT 19791. In such systems, the 
flux from small-scale sub-structures in the lensed source is se¬ 
lectively magnified by the effect of microlensing by individual 
stars in the lensing galaxy. The microlensing magnification fac¬ 
tor scales inversely proportional to the size of the source so that 
it affects mostly the flux from the smallest structures inside the 
source. Study of time and energy / wavelength dependence of 
the microlensing effect enables the study of the details of the 
source structure on the distance scales R <s: 10'® cm (or micro- 
arcseconds), comparable to the size of the supermassive black 
hole powering the AGN activity ([Torres et al.||2003) [Kochanek[ 
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[2004[ l. The microlensing data in the optical and X-ray continum 
and line emission allow to probe the structure of the accretion 
flow onto the black hole ( Eigenbrod et al.|2008[ Dai et al.|2010 


[Blackburne et al.[201 1 [[Chartas et al.[2012[ ), locate the line emis 
sion regions within the AGN central engine ( [Abajas et al. 12002} 
Guerras et al. 201 3|l and/or the probe the properties of the lensing 


galaxy ( [Mediavilla et al.|2015[ l. 

Only two strongly lensed blazars are detected in the GeV- 
TeV energy band by Eermi Large Area Telescope (EAT) ( [Abdo 
[et al.[[2015[ Cheung et al.[[2014[ l and by the MAGIC telescope 
( [Mirzoyan|2014[ ). The microlensing effect is observed in one of 
the two sources, PKS 1830-211 ( Neronov et al.|20l3] l, limiting 
the source size to be not larger than 10-100 Schwarzschild radii 
of the supermassive black hole powering the source. 

The other y-ray loud biazar which is a part of the strong grav¬ 
itational lens is B0218+357 ( [Cheung et al.[[2014| ). The system 
includes the lensing galaxy at the redshift z ~ 0.68 (Browne 
1993|l and a more distant biazar at z « 0.94 ([Cohen et al. 


et al. 


2003 . The gravitational time delay between the two images of 


the source is determined in the radio, = 10.5 + 0.4 d (Biggs 
et al. 199^, Trad = 10.1 ± 1.6 d ( [Cohen et al.[200^[Eulaers & Ma 
gain|201 n and in the y-ray band - Ty ^ 11.46 + 0.16 days (Che 


ung et al.|2014|l. The magnification factor ratio between the two 


images changes from pradio - 2 at 1.65 GHz frequency up to 
Pradio - 4 at 15 GHz, presumably due to the free-free absorption 
in a giant molecular cloud (GMC) located in the lensing galaxy 
in front of the image B ( Mittal et al.|2007| l. 

In the y-ray band, the two images of the source, separated 
by only 0.3", are not resolved. The magnification factor ratio 
is measured based on the fitting of the y-ray lightcurve of the 
source, which is composed of the contributions from the two 
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sources. This gives the result Py ~ I (Cheung et al. 2014 1 , which 
is different from the free-free absorption corrected magnification 
factor ratio in the radio band ( [Mittal et al.|2007| ). 

The difference in the magnification factor ratios between the 
radio and y-ray bands could be a signature of the microlens- 
ing ( Torres et al.||2003 Abdo et al. |2015| Cheung et al.||2014 


Neronov et al.|2015| l. This would be interesting for constraining 
the size of the y-ray source. However, the complicated chromatic 
behaviour of the radio magnification factor makes a straightfor¬ 
ward interpretation of the magnification factor difference impos¬ 
sible. 

In what follows we notice that the microlensing does affect 
the magnification factor ratio in the y-ray band. We deduce this 
from the time variability of the y-ray magnification factor ra¬ 
tio ( jCheung et al.||20T4{ [Buson et al.||201^ . Detection of the 
microlensing of the y-ray flux of the source allows us to de¬ 
rive a constraint on the size of the y-ray source. Similarly to 
the case of PKS 1830-211 ( [Neronov et al.|20r5| l, the source size 
is found to be compact, with an order-of-magnitude size esti¬ 
mate ~ lO^'* cm. The compact y-ray source size indicates that 
the y-ray emission is produced at the base of the jet, close to the 
supermassive black hole. 


2. Fermi LAT data analysis 

In this work we use the publicly available Fermi/LAT Pass 7 
reprocessed photon data se|^ The data were analysed using 
the Fermi Science Tools packag^v9r33p0, using the “Source” 
(P7REP_SOURCE) event class. The photons were selected from 
the 20° region around the position of B0218-(-357. The fluxes of 
all sources in the selected region were estimated from the likeli¬ 
hood fit, which included all the sources from the 2FGL ( [Nolan[ 
[et al.|2012) l catalogue within 28 degrees from B0218-1-357. 

The likelihood fitting procedure was repeated for each of the 
time bins of the source light curve. We used linear time binning 
of 0.5 day and extracted source fluxes in the 0.1-510 GeV energy 
band. We estimated the significance of detection in each time bin 
using the Test Statistics value (TS, [Mattox et al.[1996] ), obtained 
from the fit. Whenever the TS values dropped below 10 (which 
roughly corresponds to 3cr significance) the source detection was 
considered to be non-significant. 


2.1. Timing anaiysis 

Over the period of Fermi/LAT observations, B0218-I-357 has ex¬ 
perienced two pronounced flaring periods ( [Cheung et al.[[2014 
[Buson et al.|2015[ ). A zoom on these periods is shown in FigTjT 
The light curve of the source during the 2012 flare reveals 
strongly variable behaviour on different time scales, from the in¬ 
traday to ~ 100 d ( [Cheung et al.|2014| ). The gravitational delay 
time scale is in between these minimal / maximal scales, so that 
a straightforward decomposition of the source flux onto the two 
separate contributions from the two images is not possible. Still, 
the gravitational time delay could be reliably estimated using the 
fast variability episodes occurring on the time scale much shorter 
than Ty ([Cheung et al.|2014|l. 

The gravitationally delayed flare expected 11.5 d after the 
2014 episode is not detected. This is also clear from a detailed 
pattern of photon arrival times in the energy band above 1 GeV, 
shown in Fig. One could see that 24 photons are detected 
within a short 15 hr period of the main flare. The rise time of 


* http://fermi.gsfc.nasa.gov/cgi-bin/ssc/LAT/LATDataQuery.cgi 
^ http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/ 
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Fig. 1. Light curves of the B02184-357 in 0.1-510 GeV band over its 
2012 (top) and 2014 (bottom) flaring episodes. The bins size is 0.5 days 
for the 2012 and 1 day for 2014 light curve. Blue and grey data point 
mark the significant and non-significant detections correspondingly. 
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Fig. 2. Arrival times and energies of photons above 1 GeV during the 
2014 flare (bottom panel) and during the period after the gravitational 
time delay. Grey vertical bands mark the Good Time Intervals of Fermi 
/LAT. 


the flare is comparable to the time interval between subsequent 
exposures of the source by the LAT, ^ 3 hr. Detection of the 
delayed flare with comparable flux and comparably sharp rising 
front would provide a precision measurement of the gravitational 
time delay with the minimal possible error bar of 3 hr. However, 
no delayed flare is detected as it is clear from the upper panel 
of Lig. ^ down to the flux level which is an order of magnitude 
lower than the flux of the main flare. 


2.2. Gamma ray magnification factor ratio 

Lstimation of the magnification factor ratio between the two im¬ 
ages of the source is difficult in the y-ray band because the lensed 
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images of the source are not resolved. As a result, the light curve 
of the y-ray source has contributions from both images; 


F,ot^ ^^yF{t)+F{t-Ty) ( 1 ) 

Here denotes the magnification ratio of the leading image and 
Ty is the gravitational time delay. 

In real y-ray data analysis, the continuous function F{t) is 
replaced with a set of measurements of the source flux at discrete 
moments of time Fi - F(ti) and Eq.[T]becomes: 

Ff = ^iyFi + (2) 


where A denotes the shift Ty in terms of the light curve bins size. 
This equation can be re-written in the matrix form: 
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The source flux at the start of the observations {tstart) is com¬ 
posed of the total flux of the two images in the time interval 
preceding the tstart- In this way the system in Eq. can only 
be solved if the boundary condition for the moments of time 
t < tstart is provided. In the case of isolated flares one can as¬ 
sume that the source flux before the moment tstart is equal to a 
constant quiescent flux Fbase, as it is done in Eq. 0. 

If Ty and fiy are known, an approximate solution of Eq 
be found from the least-square minimization under the condTtion 
Fi > 0. The uncertainties of the flux estimates can be propa¬ 
gated to this system by simply dividing each value of E'"' and 
/-th row of the matrix by the corresponding uncertainty cr,. This 
transforms the least-square solution of the system in Eq. |^to the 
least-;^^^ one. 

The value of the time delay Ty can be found from the com¬ 
bined source lightcurve using e.g. the autocorrelation function 
analysis, which results in Ty a; 11.5 days (|Cheung et al.|2014|l. 
The magnification factor ratio fiy is more dithcult to hnd because 
of the large number of uknowns E, in the system Q, equal to the 
number of equations. To assist the situation, an additional con¬ 
straint on the solution F, can be imposed. 

In the absence of accidental fine tuning of parameters, the 
true solution E(f) of Eq. (or Eq. should not contain an au¬ 
tocorrelation at the time scale of the gravitational delay. In this 
way the best solution of the system in Eq. |^can be found by 
minimizing the residuals between the autocorrelation value of 
E(f) around the position of expected time delay Ty and its local 
approximation, e.g. by a power law. 

In a particular case of B0218-1-357, the major source flare 
in 2012 was preceded by the period of a relatively constant 
flux (Cheung et al. 20141 1, which allows the estimation of the 
baseline flux Fbase, needed to solve the system 0- 

Solving the system (|^ using the procedure described above 
and assuming fiy constant in time we And that no acceptable 
solution could be found. To get an acceptable solution, the as¬ 
sumption on the constancy of jUy has to be relaxed. In principle, 
the time scale of variability of fiy is not known a-priori. How¬ 
ever, general considerations of microlensing suggest that even if 
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Fig. 3. Top left: Decomposition of the different flaring episodes of the 
B0218-)-357 major flare in 2012. Top right: Derived magnification factor 
ratio for each of the flaring episodes. Black dashed line marks the ap¬ 
proximate value of the magnification factor ratio in radio band prad ~ 4. 
The red line depicts the tentative association of the observed magnifica¬ 
tion factor behaviour with one the caustics crossing events in our sim¬ 
ulations, computed for Rsource ~ 10*'* cm. This event is shown in lower 
panels. Bottom: Microlensing caustics magnification patterns, that can 
be associated with the detected variability of magnification factor over 
the 2012 and 2014 flares. Distances are given in the source plane. The 
patters are shown separately for the leading (left) and delayed images 
(right). Solid white lines mark the putative trajectories of the images 
during the 2012 flare, dashed lines depict the displacement of the im¬ 
ages over the 2012-2014 period. 


the source size is as small as the size of the supermassive black 
hole, the variability time scale of the magnification factor ratio 
is the months-to-years range. 

The 2012 major flare could be split on several shorter flar¬ 
ing episodes, identified by [Cheung et al.| ( |2014) l. To study the 
variability of py, we modify our procedure for solving the sys¬ 
tem ([^ by allowing py to change from sub-flare to sub-flare (in 
a step-wise manner). We search for a solution of the system ([^ 
performing a scan over the magnification factor ratios for each 
flaring episode. This results in an acceptable solution E(f) shown 
inEig.|^ 

To estimate the uncertainties on the derived best-fit py for 
each flaring episode, as well as on the intrinsic lightcurve E(f), 
we have simulated a number artificial light curves E,o,(f) in 
which the flux values in each time bin are randomly scattered 
within the error bars of the original flux measurement. These 
light curves were substituted to the system]^ which was solved 
as described above. The distributions of the best-fit py for each 
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Fig. 4. Distributions of the best-fit values of the magnification factor 
ratio p for the three brightest flaring episodes in 2012 (3 central panels 
in Fig.[^, derived through the Monte Carlo simulations. 


flaring episode is shown in Fig.|^ In order to quantify the scatter 
of the derived values of the magnification factor, we have fitted 
these distributions with a Gaussian profile. The width of the pro¬ 
files is shown as the errorbars on the measurements of py in each 
episode in the right panel of Fig.|^ 

We would like to note, that a similar behaviour of the mag¬ 
nification factor during the 2012 flare can be also seen in the 
analysis of |Cheung et al. ( |2014| l, although the superposition of 
the leading and delayed flaring components did not allow there 
to make firm conclusions. 

Contrary to Cheung et al. ( |2014| l, we do not find an evidence 
for the “orphan” sub-flares with no associated delayed emission. 
Detection of such sub-flares would imply a very rapid variation 
of the magnification factor, which would be interesting for the 
search of the microlensing. However, we find that the data are 
consistent with relatively slow modulation of the magnification 
factor ratio over the entire 2012 flare. 

The magnification factor ratio, measured during the second 
flare of the source in 2014, shown in Fig.[T] can not be directly 
established, as the source flux at the moment Tfiare+Ty is consis¬ 
tent with the flux in-between the putative flares, i.e. the delayed 
emission was not detected by Fermi/LAT. The rising part of the 
flare in the energy band above 1 GeV is very sharp, with the 
overall duration of several hours, see Fig. During the 15 hr 
flaring time, 24 photons were detected within the circle of the 
radius 0.5° around the source. At the same time, only two or 
three photons were detected during the 15 hr time window of 
the delayed flare (compare the lower and upper panels of Fig.|^. 
It is not clear to which of the two images belong the two / three 
photons detected during the delayed flare time window. The ratio 
of the photon counts of the “prompt” and “delayed” 2014 flares 
provides a lower bound on the magnification factor ratio Py > 5 
during this flare. 


3. Discussion 


The analysis of the previous section shows that the magnifica¬ 
tion factor ratio py has varied during the 2012 flaring period and 
has further changed from 2012 to 2014. An immediate conse¬ 
quence of the detection of time variability of py is the conclu¬ 
sion on the detection of the microlensing effect in the y-ray band. 
An alternative possibility for the variations of the magnification 
factor ratio (as discussed by Bamacka et al. ( 2015| l in the case 
of PKS 1830-211) would be different locations of the 2012 and 
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2014 flaring regions inside the source. This, however, would re¬ 
quire the displacements of the order of the Einstein radius of the 
(macro)lens, ~ 1 kpc, even during one and the same flaring pe¬ 
riod. The scale of displacements is much larger than the typical 
estimates of the sizes of the y-ray emission regions in blazars 
/ quasars and radio galaxies and, moreover, in the case of the 
2012 flaring episode, it is looks unlikely that two distant flaring 
regions would flare simultaneously within a 100 d period. 

On the other hand, the assumption of microlensing provides 
a natural explanation to the detected magnification factor ra¬ 
tio behaviour. Using the simulated micorlensing caustic patterns 
maps, described in more detail in section |3.2[ we were able to 
tentatively associate the variation of py over the 2012-2014 flares 
with the event of caustic crossing of the delayed image in 2012 
and a change in magnification of the leading image between 
2012 and 2014 due to the movement over a larger scale pattern 
of lensing caustics. This tentative association is illustrated in the 
lower panels of Fig. and provides a reasonable description of 
the observed behaviour. 

3.1. Qualitative estimates of the microlensing constraint on 
the source size 

The microlensing magnifies the flux of the source if it is smaller 
or comparable in size to the lens Einstein radius: 


Re(M) = 


4GM Dls 
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where Ds,Dl,Dsl are the angular diameter distances to the 
source, the lens and between the source and the lens, D is the 
overall distance scale and M is the mass of the lens. 

If the lensing centres in the lens galaxy are stars, only a 
source of the projected size less than Re corresponding to the 
lenses of the masses ~ Mq is subject to the microlensing. In par¬ 
ticular, the radio source which is known to be large (it is resolved 
in the radio band) Rradh ~ 10 pc ( Mittal et al.|2007[ ) is not influ¬ 
enced by the microlensing by stars. 

The fact that the micorlensing effect is detected in the y-ray 
band immediately implies a constraint on the size of the y-ray 
source projected on the plane of the sky: 


Ry.proj < Re{Mq) 3 X 10 


16 


(5) 


under the “default” assumption M ~ Mq. 

The microlensing magnification factor is variable in time due 
to the relative motion of the source, the lensing centres and the 
observer. Maximal magnification is achieved at the moments of 
crossing of the microlensing caustics in front of one of the two 
images of the source. At these moments, an order-of-magnitude 
estimate for the magnification factor is (Chang|1984|l: 


- ^RE/Ry.proj - 


10 




3 X 10^"^ cm 


(6) 


Assuming that the high-frequency (15 GHz) radio observations 
provide a fair estimate of the macro-lensing magnification fac¬ 
tor ratio, pradio - 3.7, the overall magnification factor ratio is 
expected to vary in the range Pradiolp"""° < Ry < PradioP"""”, 
where we have made a simplifying assumption that both images 
of the source are affected by the microlensing in a similar way. 

Measurement of the variability of the magnification factor 
ratio from py < 0.4 + 0.1 to //y > 5 allows to make an estimate 
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of jUmicro ~ 10. This gives an oder-of-magnitude estimate of the 
source size 
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y,pro] 


3 X 10 
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10 


cm 


(7) 


An additional constraint on the projected size of the source 
stems from the measurement of the time scale of the variabil¬ 
ity of /ty. The variability shown in the right panel of Fig. [^is 
characteristic for a caustic crossing episode during which 
reaches a maximum / minimum. The typical timescale of the 
caustic crossing events is determined by the relative velocity v 
of the source with respect to the caustics pattern and the ob¬ 
server ( |Kayser et al.|198^ . Assuming that the velocities of the 
source and the lens are of the same order as the velocity of the 
Sun w.r.t. the Cosmic Microwave Background reference frame 
V ^ 300 km/s ( |Kogut et al.|1993| ) one could find that the strong 
magnification events related to the crossing of the source by the 
microlensing caustics last for 


Aty — Ry,pro j rel ~ 100 




3 X 10'"^ cm 


Vrel 


300 km/s 


( 8 ) 


which is consistent with the data of the 2012 flaring episode. 

The caustic crossing also provides a natural explanation for 
the exceptional brightness of the 2012 flaring episode. In the 
absence of the microlensing, the flux from the dimmer image 
of the hlazar is by a factor of three lower than the flux of the 
brighter component. A moderate flaring activity typically pro¬ 
duces the flux at the level F ~ 10“® ph/(cm^ s) in the brighter 
image of the source and just at the level F ~ 3 x 10“^ ph/(cm^ 
s) in the dimmer image. An order of magnitude change in the 
magnification factor of the dimmer image, caused by the caus¬ 
tic crossing resulted in a boost of the dimmer image flux up to 
F ~ 3x 10“® ph/(cm^ s), so that the overall source flux has grown 
up to F tot ~ 4 X 10“^ ph/(cm^ s) during the flaring period. 


3.2. Numerical modelling of the microlensing magnification 
factor statistics 

Qualitative estimates presented above are uncertain because of 
the low statistics of the microlening events (caustic crossings) 
in the Fermi / LAT data. Only one such event is detected after 
six years of Fermi / LAT exposure. The maximal magnification 
factor and the time scale of the caustic crossing events fluctuate 
due to the complex shapes of the caustic patters. These fluctu¬ 
ations introduce an uncertainty into the estimates of the source 
size based on the microlensing data. 

This uncertainty is illustrated by Fig. where the expected 
statistics of the magnification factor ratio with account of the 
microlensing effect is calculated for different y-ray source sizes. 
The calculation presented in this figure was performed via the 
simulation of the microlensing magnification maps using the 
inverse ray shooting method ( |Kayser et al. |1986[ [Schneider &| 
|Weiss|1986||1987| l. The simulated region of space encompasses 
0.2x0.2 pc in the source plane and contains 300-1700 stars- 
lenses, depending on the assumed optical depth with respect 
to microlensing, which was scanned in range Z = [0.2,1.0]. 
The lenses mass distribution was chosen to follow the IChabrier] 
( [2003 | l initial mass function with the cut-off at 10 Mq. For each 
simulation run we have randomly selected 10 trajectories for 
each of the two lensed images of the source. The “magnifica¬ 
tion light curves”, computed along these trajectories, represent 




Fig. 5. Top: Distribution of the average durations of the microlensing 
magnification episodes obtained from the simulations for E = 0.5 and 
V = 300 km/s. Horizontal shaded region depicts the range 60-120 days 
of the apparent variability time scales of the magnification factor ratio in 
B0218-)-357 and the vertical one depicts the microlensing magnification 
range, suggested by the observations. Bottom: Cumulative probability 
of observing microlensing magnification > p for different source sizes, 
estimated from the simulations. The vertical shaded region is the same 
as in the upper panel. 


the magnification of each of the images during the source move¬ 
ment with respect to the lens. In order to compare with the de¬ 
tected magnification factor ratio, we have computed the ratio of 
the trajectories for both images. Using these simulated magnifi¬ 
cation ratio curves, we then calculated the average durations of 
the episodes of high magnification > fi. 

Fig.|5] shows the dependence of the variability timescale of 
the microlensing on the source size, obtained from the sim¬ 
ulations. One could see that a source of the size larger than 
3 X 10'^ cm is typically not expected to change the magnification 
factor ratio on the time scale shorter than two years, as observed 
in B0218-1-357 system. At the same time, a source of the size 
3 X 10*^ cm would exhibit variations of the magnification factor 
ratio already on the time scale of the 2012 flare. Thus, the detec¬ 
tion of the microlensing implies an estimate of the source size 
Ry<3x 10'4 (v/300 km/s) cm. 

The largest uncertainty of the estimate stems from the un¬ 
known projected velocity of the source. In principle, this veloc¬ 
ity could be up to V ~ c, because the flaring source could be 
a relativistically moving blob inside the jet. Such fast motion 
of the flaring blob would explain the “orphan” high magnifica- 
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tion sub-flares tentatively found in the analysis of |Cheung et al.| 
( 2014| l. At the same time, this would imply a significantly larger 
source size ^ 3 x 10*^(v/c)(Afy/100 days) cm in order to 
account for the detected duration Afy ~ 100 days of the high 
magnification event during the 2012 flaring episode. The possi¬ 
ble apparent superluminal motion v > c of the highly relativistic 
blobs would make this estimate even larger. Such source size is 
much larger than the size of the Einstein ring of the micro-lens 
in B0218-I-357. In this situation, the microlensing would not be 
able to significantly affect the source flux, as can be seen from 
the upper panel of Fig. Still, a mildly relativistic source with 
V ~ 0.1c of the size < 10'® cm would already occasionally pro¬ 
duce high-magnification events with ~ 10. 

The cases of large fast moving source and smaller slow mov¬ 
ing source could be distinguished based on the expected statis¬ 
tics of the high-magnification events. During more than six years 
of Fermi / FAT observations, only one caustic crossing / high 
magnification event is detected. This fact provides a rough esti¬ 
mate of the probability to find the value > 10 in the mi¬ 
crolensing lightcurve, > 10) ^ (60 d/ (6 ■ 365 d) ^ 3%. 

Fower panel of Fig. shows that this estimate favours the com¬ 
pact source of the size < 3 x 10'"' cm. If the source size is larger, 
the fraction of caustic crossing events with sufficiently high mag¬ 
nification decreases. Of course, with only one detected caustic 
crossing event, the estimate of > 10) ^ 3% is not very 

precise. Further monitoring of the source is needed to confirm 
the hypothesis of the small source size. 


3.3. Comparison of the microlensing constraint with the 
constraint from the minimal variability time scale 

The microlensing constraint on the source size is consistent with 
the constraint stemming from the minimal variability time scale. 
The 2014 flare was particularly sharp in time, with the overall 
duration Tfiare ^ 15 hr in the energy band above 1 GeV (see Fig. 
1^. The rise time of the flare is comparable to the interval be¬ 
tween the subsequent Fermi/FAT exposures of the source spaced 
by the Trise - Tlat - 3.2 hr period. This fast variability con¬ 
strains the linear source size to be (|Celotti et al.|1998[ |Neronov| 
let alJ20081 l 


Ry < 3 X IO'^'F^ 


^rise 
104 S. 


cm 


(9) 


where F is the bulk Forentz factor of the jet. If the jet is aligned 
along the line of sight within a viewing angle 0 ~ F^', the size of 
the emission region projected on the plane of the sky is limited 
to 


Ry.proj = Ry6 < CT„i J 3 X 10 


F 

rise 

10 

[l04 sJ 


( 10 ) 


Assuming the jet bulk Forentz factor F ~ 1-10, one could 
find that the expected projected size of the y-ray emission re¬ 
gion derived from the variability analysis is comparable with the 
estimate of the size of the y-ray emission region from the mi¬ 
crolensing analysis. 


pairs in interactions with UV photons of energies about e ~ 
10 eV and, thus, might be significantly absorbed if the UV pho¬ 
ton density is sufficiently big. 

Fumin osity of the source in this energy range is about Luv ~ 
1044 erg/s (Abdo et al. 2010| . If this luminosity stems from the 
accretion flow, the UV photon distribution is isotropic in the 
black hole reference frame. If the UV source size is Ruv^ the 
number density of the UV photons is about 


^ph 


~ 10" 

Ruv 
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47TRlyCe 

10'® cm 


1044 erg/s 
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The mean free path of the 100 GeV y-rays is limited to 

1 ... KjI\/ “ I ^TTf/ 
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( 12 ) 


and the optical depth of the source with respect to the pair pro¬ 
duction is 


iyy 



Ruv 

-1 

Ruv 

Ayy 

o 

o 
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1044 erg/s 


(13) 


Thus, if the UV source is as compact as the y-ray source, R/jy < 
lO'® cm, the highest energy (100 GeV) photons could not escape 
from the source. Detection of such photons by MAGIC telescope 
(|Mirzoyan|2014)l imply that the UV source is large. 


Ruv ^ 10 


18^-1 

rr 


Ruv 

1044erg/s 


cm. 


(14) 


The origin of the source flux in the UV band is uncertain. 
The source spectrum does not exhibit a UV bump characteristic 
for the black hole accretion disk. If the emission is dominated 
by the jet, the bulk of the UV flux is most probably produced 
via inverse Compton scattering by relatively low energy elec¬ 
trons, so that the large size of the UV source does not appear 
surprising. Otherwise, if the accession flow onto the black hole 
in B0218-1-357 is radiatively inefficient, the UV emission is pro¬ 
duced by the inverse Compton scattering by the non-relativistic 
electrons at large distances from the black hole. In this case the 
large size of the source is also naturally expected. 

The constraint on the UV source size could be directly veri¬ 
fied via the (non)-observation of the effect of the microlensing in 
the UV band. Indeed, a source of the size Rjjy > Re(Mq) should 
be almost not affected by the microlensing. Simultaneous y-ray 
and UV observations of the source could be used to measure the 
time variability of the magnification factor in the UV band (or the 
absence of it). Variations of the magnification factor ratio com¬ 
parable to those observed in the y-ray band should produce an 
order-of-magnitude variations of the UV flux. Tracing the varia¬ 
tions of the magnification factor ratio via y-ray observations and 
comparing them with the variations of the UV flux should pro¬ 
vide a test of the presence / absence of the microlensing effect in 
the UV band. 


4. Conclusions 


3.4. Problem of escape of very-high-energy j-rays from the 
source 


Compactness of the y-ray emission region in B0218H-357 might 
pose a problem of escape of the very-high-energy y-rays, which 


have been recently detected by MAGIC telescope (Mirzoyan 
2014 1 . Gamma rays with energies above 100 GeV produce 


We have demonstrated that the y-ray signal from the gravitation¬ 
ally lensed blazar B0218-I-357 is affected by the microlensing. 
Our conclusion is based on the detection of variability of the 
magnification factor ratio in the range py ^ 0.4 to py > 5 both 
on short (< 100 days) and long (~ 600 days) time scales. 

Using the data on the time variability of the magnifica¬ 
tion factor ratio, we have derived a constraint on the projected 
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size and location of the y-ray emission region. The order-of- 
magnitude estimate of the size is Ry ~ 10'^ - lO'^ cm. The 
microlensing estimate of the source size is consistent with an es¬ 
timate stemming from the minimal variability time scale of the 
source. The y-ray emission region is not moving with relativistic 
speed so that it’s most probable location is close to the hlazar 
central engine, at the base of the jet. 

Constraints on the size and location of the y-ray emission 
region in B0218H-357 are similar to those derived from the 
detection of microlensing in the other strongly lensed hlazar 
PKS 1830-211 ( [Neronov et al. 12015) 1. This demonstrates that the 
small size and location close to the supermassive black hole are 
generic features of y-ray emission from blazars. 
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